Slow quench dynamics of a trapped one-dimensional Bose gas confined to an optical lattice 
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We analyze the effect of a linear time- variation of the interaction strength on a trapped one-dimensional Bose 
gas confined to an optical lattice. The evolution of different observables such as the experimentally accessible 
onsite particle distribution are studied as a function of the ramp time using time-dependent exact diagonalization 
and density-matrix renormalization group techniques. We find that the dynamics of a trapped system typically 
display two regimes: for long ramp times, the dynamics are governed by density redistribution, while at short 
ramp times, local dynamics dominate as the evolution is identical to that of an homogeneous system. In the 
homogeneous limit, we also discuss the non-trivial scaling of the energy absorbed with the ramp time. 
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Manipulating many-body quantum systems by time- 
varying their control parameters is a practical challenge of 
technological importance in many areas of physics including 
condensed matter, quantum information, and cold atomic and 
molecular gases. However, our understanding of the quantum 
dynamics of many-particle systems and the identification of 
their universal dynamical features is still in its infancy [1, 2]. 
In recent years, it was suggested that the Kibble-Zurek mecha- 
nism J^t], originally developed to describe the evolution of the 
early universe, could explain the dynamics of systems across 
quantum phase transitions. Despite a few successes, the va- 
lidity of this theory to describe the evolution of all quantum 
systems is still not accepted. Unbiased theoretical methods, 
going beyond scaling arguments, are required to understand 
the dynamics of both homogeneous and inhomogeneous non- 
integrable quantum systems. In an attempt to shed some light 
on the evolution of bosonic systems subjected to a change of 
their control parameters, experiments on 4 He [4] and more 
recently on cold atoms were reported. Quenches, conducted 
on a trapped bosonic quantum gas loaded into an optical lat- 
tice, were performed by changing the depth of the lattice over 
a given time interval j5Ha]. In these experiments, the out- 
of-equilibrium processes were investigated by considering the 
behavior of local observables such as the density, compress- 
ibility and onsite particle distribution. For slow to moderate 
quenches, two different evolution regimes which depend on 
the ramp time and on the initial interaction strength were ob- 
served. These two regimes are believed to be related to the 
local and global dynamics of the system How well 

these experiments can be used to clarify the universal dynam- 
ics of homogeneous systems has not been addressed yet. 

In this work, we provide answers to this question by analyz- 
ing the response of bosons stored in a one-dimensional (ID) 
lattice to a slow increase of their interaction strength using the 
unbiased methods of exact diagonalization (ED) and density- 
matrix renormalization group (DMRG). In the homogeneous 
case, we find, in addition to the sudden quench and quasi- 
adiabatic behaviors observed for fast and slow ramp times re- 
spectively, that at intermediate ramp times the absorbed en- 
ergy scales non-trivially with the ramp duration. In the pres- 



ence of a trap, we identify two distinct regimes as a function of 
the ramp time. For long ramp times, the evolution is governed 
by density redistribution, whereas for shorter ramp times, the 
evolution is dominated by intrinsic local dynamics and mass 
transport is absent. This last response is the same as that of an 
homogeneous system. 

We carry out our study using the ID Bose-Hubbard model: 
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with bj the operator creating a boson at site j, nj = bjbj 
the density operator and J and U the hopping and interaction 
amplitudes. The chemical potentials account for an exter- 
nal confinement. At commensurate fillings, a quantum phase 
transition from a gapless superfluid (SF) to a gapped Mott in- 
sulating (MI) state occurs (at U w 3.3J for n = 1 iflfjh . 
The slow quench is performed by increasing the interaction 
strength linearly, i.e. U(t) — Ui + ^SU with r the ramp 
time, SU = U '/ — Ui the quench amplitude, and ?/»(/) the 
initial (final) interaction. This can be achieved experimen- 
tally using a suitable Feshbach resonance 111 111 . Aspects of 
linear quenches have been discussed previously using various 
approximate methods . Here, time evolution is com- 

puted numerically on chains of size L, using both ED with pe- 
riodic boundary conditions and an onsite boson cutoff M > 7, 
and DMRG [13] with open boundary conditions and M = 6. 
The convergence of the DMRG results both with the num- 
ber of states (a few hundreds) of the reduced space and the 
time-step of the Trotter- Suzuki time-evolution were checked. 
Denoting by E oi /f the initial and final ground-state energies 
and Ef the energy obtained at time r, we introduce the heat 
as the energy absorbed by the system: Q = Ef — Eqj. Note 
that we only consider the time-evolution during the ramp and 
not the relaxation once the ramp is completed. Further, the 
derivative of the chosen ramp has a discontinuity at the initial 
time, by which higher modes might be excited. 

Homogeneous system - We first aim here at understanding 
the intrinsic evolution of local observables by studying the ho- 
mogeneous limit as many features identified in such systems 
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with the characteristic ramp time r c given by 
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FIG. 1. (Color online) Ui = 2J and £// = 3 J at n = 1. (a) Scaling 
of the heat put into the homogeneous system vs inverse ramp time r 
showing a non-trivial exponent in the near-adiabatic regime. Inset : 
collapse of data for different SU (ED), (b) Comparison to perturba- 
tion theory (r c ~ 3.47ft/ J). Inset: Spectral function associated with 
perturbation theory for U = 2 J (ED). 



are also relevant to the trapped case. In Fig. [Tf a X we show 
the heat produced by the quench for a given SU as a function 
of the inverse ramp time r _1 . In finite-size systems, three 
regimes are typically observed: at large-r (adiabatic limit), 
Q ~ r~ 2 with oscillations [1], here associated with a finite- 



size gap Hl4jl : at intermediate r, we find a non-trivial power- 
law behavior; finally, for short r, Q approaches the sudden 
quench limit II 1511 quadratically. 

These results can be compared with time-dependent pertur- 
bation theory in the initial Hamiltonian basis {\n)} where the 
"ramp velocity" v = SU /r is the small parameter. To second 
order in v, the energy variation reads: 



AE(t) = vtB 
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with B = J2j n ji n j ~ the perturbation operator and 
B n = (n\B\0). Numerically, we access the \B n \ 2 and ener- 
gies oj n — E n — E in the low-energy regime using the Lanc- 
zos algorithm. Considering only the final time r and introduc- 
ing the energy spectral function, B(lo) — J2 n |-B„| 2 <5(u;— cj„), 
one gets for the heat 

Q(t,SU) = Q°(SU) 
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E j is the exact sudden 
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where Q°{SU) = SUB + E ,i 
quench expectation. Counter-intuitively, even though this re- 
sult is perturbative in v, it yields the correct result in the large 
velocity v limit as the ramp time r becomes short enough in 
this regime. Thus, this perturbation theory provides short-r 
corrections away from the sudden quench limit, which can be 
calculated using ground-state observables. Indeed, assuming 
that B(lj) has a support with an upper bound or decays expo- 
nentially, ones finds that Q(r, SU) = Q°(SU)-L^{t/t c ) 2 
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where K is the kinetic term ■ +h.c.]. The magnitude 

of t c is consequently directly connected to the equilibrium 
three-point correlator B[B, K\. The scaling of t c with L de- 
pends on the typical behavior of the correlator with distance. 
In our model, we find with negligible finite-size effects that 
t c is an increasing function of Ui, indicating that the correla- 
tor drops off rapidly for any U. Using Eq. (Q]), one obtains a 
quantitative agreement up to intermediate velocities and close 
to the power-law regime (see Fig.QJb)). Furthermore, the de- 
tails of B(lv) do not alter much this short-r regime as a trun- 
cated triangular approximation of B(lo) ("model" in Fig.QJb)) 
reproduces well this regime. 

We observe that the accuracy of Eq. (Q~|) improves with de- 
creasing SU. In the small SU limit, the quench essentially 
becomes a probe of the initial ground-state dynamics. To sec- 
ond order in SU, the maximum heat behaves as Q°(SU) ~ 
SU 2 2j„>o \B n \ 2 /u n . Combining this result with Eq. (Q]), we 
find that Q(t,SU)/Q°(SU) ~ /(t) is a function of r only. 
In the Inset of Fig. \V[a), we see that curves do collapse well 
onto each other for SU < 0.6; while they do not for larger 
SU as in this case heat production depends on both r and SU. 
Remarkably, the small SU limit also provides information on 
the large-r (adiabatic) regimes as, in that case, the first part 
of the integral in Eq. (fTJ cancels with Q°(SU). There, we 
recover the results of Ref. 16. In particular, we see that for 
a gapped spectrum Eq. (fl} reproduces the typical r~ 2 decay 
combined with some oscillating terms. Lastly, we discuss the 
exponent rj w 1.35 observed for Ui = 2 J and SU = J. For 
large system sizes, this behavior is found above the sudden 
quench limit but below a ramp time set by the inverse finite- 
size gap fl7ll . We cannot relate this exponent to various pre- 
dictions on approximated versions of the model fljj] , nor to 
the above perturbative approach. Similar non-universal be- 
haviors were recently reported in Ref.[l7t Moreover, the Inset 
of Fig. Q2a) shows that this exponent depends on SU. This 
may have different reasons. First, contributions of intermedi- 
ate interaction values at which the system possibly possesses 
different low-energy physics might play a role. Second, a dis- 
continuity in the derivative of the considered ramp form might 
cause higher energy excitations. Last, energy redistribution 
between excited states, which we found to be negligible only 
for small quenches, could also play a role. This exemplifies 
that, typically, only the small SU limit can display universal- 
ity. 

Following experiments, we turn to the discussion of local 
observables. kj = (n 2 ) — (rij) 2 denotes the compressibility 
at site j while P n is the probability of having n bosons onsite. 
The evolution of Fig.|2ja) starting from the SF phase typically 
displays the same three regimes as for the heat, with finite- 
size oscillations close to the adiabatic limit. Entering the MI 
at n = 1, Pi increases with U while other P n concomitantly 



FIG. 2. (Color online) Observables for a quench from Ui = 2J to 
Uf — 4 J (ED) (a) Occupation probability P„ and compressibility k 
vs ramp time r in the homogeneous system, (b) Perturbation theory 
of Eq. ([2} gives a quantitative prediction for L — 11. The L = 16 
data gives an idea of finite-size effects. 

decrease. The perturbative prediction for a real symmetric ob- 
servable A reads, at time r and to first order in v, 

A(r, 6U)=A -2 6 -^Y: ^ " t KT) ^ ■ (2) 

n^O ™ 

Here again, it is remarkable that the simple perturbation the- 
ory provides a quantitatively good (up to a percent in this case) 
prediction of the exact time-evolution (see Fig. [2b)). 

Trapped system - In current cold atom experiments, an ex- 
ternal parabolic trap is usually present, i.e. fij = —V(j — 
-^ii) 2 . As the density increases from the edges to the center, 
for large U, MI domains coexist with SF regions 11811 . Thanks 
to recent advances, measurements with single-site resolution 
are now possible in 2D cold atoms lattice setups 

0,11, 

en- 
abling one to focus on a particular domain. However, in an 
out-of-equilibrium situation, density gradients and parameter 
changes induce flows of particles. Consequently, the density 
redistribution will have an impact on locally measured quan- 
tities. These transport phenomena occur on timescales that 
depend on the velocity of excitations (typically controlled by 
J) and the domain sizes. 

Generally, the timescales for density redistribution are 
longer than the intrinsic dynamics of a local observable. We 
show in this section that, in the fast quench regime, an in- 
trinsic evolution of the observables occurs while the density 
redistribution remains frozen. On the contrary, this redistri- 
bution dominates in the slow ramp regime. To demonstrate 
this point, we systematically compare the evolution of observ- 
ables in the trapped cloud with its homogeneous counterpart 
(choosing the same initial density in the center). We consider 
the evolution of a trapped cloud whose majority of atoms are 
initially in a SF state as the system features no or only weak 
MI regions signaled by a trough in the compressibility. Fig. [5] 
displays comparisons between observables taken at the cen- 
tral site as a function of the ramp time, and for two different 
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FIG. 3. (Color online). Comparing local observables at the center of 
the system between homogeneous and trapped (V — 0.006J) con- 
figurations as a function of the ramp time r (DMRG calculations). 
The grey area shows the timescales over which the responses are 
identical, (a) slow change from Ui — 2 J to Uf = 4 J. (b) Ui = 4 J 
to Uf = 6 J. 



interaction quenches: (a) U = (2 -> 4) J, (b) U = (4 — !• 6) J. 

For case (a) of Fig. [3] we see that for short ramp times the 
central density does not rearrange. In fact, the density remains 
constant for all ramp times shorter than r m h/J. Meanwhile, 
the compressibility k and probabilities P\ii show well pro- 
nounced variations. Remarkably, for this range of r, we find 
all observables to be in excellent agreement with the homo- 
geneous system predictions. This result supports our previ- 
ous statement that intrinsic local dynamics dominate at short 
ramp times. For longer ramp times t > 2hj J, we find a 
clear change in the central density, associated with the onset of 
particle currents 111911 . Naturally, this reduction of the density 
modifies the particle distribution and compressibility, driving 
them away from the constant density behavior and close to the 
adiabatic expectation for the trapped cloud (with slight oscil- 
lations). On Fig. Ufa), we show the actual time evolution of 
the density and compressibility corresponding to a ramp time 
t = 20K/ J. The density configuration remains almost frozen 
up to time t — 5H/ J even though the corresponding ground 
state considerably differs, and begins to evolve afterward. At 
later times, the density profile broadens reaching a wider size 
than the one expected for the ground state at the correspond- 
ing U/J. At the end of the evolution process, the system is 
therefore in an excited state. In comparison, the compress- 
ibility distribution evolves much faster than the density. At 
each time step, the compressibility is relatively close to its 
corresponding ground state value. The difference between the 
time-evolved and ground-state values are most likely due to 
the unrelaxed density profile. Consequently, this direct time- 
analysis of a slow quench reinforces our previous conclusion 
that the intrinsic evolution of the observables and the density 
redistribution are characterized by two different timescales. 

For case (b) of Fig. [3] the first striking feature is the even 
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FIG. 4. Slow quenches with r = 20ft/ J, L = 64 and V = 0.006 J. (a) f/ = (2 -+ 4) J, (b) U = (4 -> 6) J. Contour plots: density n 
(up) and compressibility k (down) as a function of position and time. Graphs: n (up) and k (down) vs position for various times t (full lines). 
Circles and squares stand for ground-state predictions with corresponding U (t). 



slower evolution of the central density with increasing ramp 
times when compared to the situation discussed before. Up to 
r w 2H/J, the density remains very close to its initial value 
and the other local observables are identical to their corre- 
sponding homogeneous value. For longer ramp times, the cen- 
tral density then shows an oscillatory decrease and the other 
observables deviate from their constant density counterparts. 
However, even for the longest ramp times considered, the den- 
sity has still not reached the adiabatic regime. The slowdown 
of the density rearrangement is attributed to the emergence of 
Mott "barriers" in regions where rij « 1. This effect is un- 
derlined in the snapshots of Fig. Hfb). These "barriers" arise 
due to the rapid reduction of particle fluctuations at large U 
and the local reduction of transport. On Fig. Sfb), we ob- 
serve a much slower redistribution of the density than for the 
Ui = (2 — ► 4) J evolution in (a). In this case, the density 
profile is actually narrower than the ground state profile. The 
compressibility presents the same fast evolution at short time, 
building "barriers" before the currents can establish and pre- 
cluding the redistribution from taking place. 

In summary, we performed a detailed study of slow 
quenches in the trapped ID Bose-Hubbard model. We iden- 
tified two dynamical regimes. For short ramp times during 
which currents cannot develop, the intrinsic local dynamics 
dominates and the response is equivalent to that of an homo- 
geneous system. Many features of this regime are well de- 
scribed by perturbation theory. For longer ramp times, cur- 
rents do set in and the dynamics is governed by non-trivial 
transport phenomena. For large final interaction strengths, the 
global timescales are significantly enhanced by the presence 
of "Mott barriers". Our results agree qualitatively well with 
recent experiments focusing on transport dynamics of cold 
atoms across the SF-MI transition |g|8|]. Indeed, based on 
our results, we can argue that the two experiments were prob- 
ably conducted in different regimes. In Ref.Q the slowdown 
of dynamics, due to the presence of Mott regions, indicates 
that their experiment was carried out in the regime with global 
dynamics. In contrast, the short timescales observed in Ref. 8 



are most likely related to local dynamics. 

Note added: During the final stages of the preparation of 
this manuscript, a mean-field study of the 2D version of the 
model appeared on arXiv ll20ll . 
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